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We focus on determining the separability of an unknown bipartite quantum state p by invoking 
a sufficiently large subset of all possible entanglement witnesses given the expected value of each 
element of a set of mutually orthogonal observables. We review the concept of an entanglement 
witness from the geometrical point of view and use this geometry to show that the set of separable 
states is not a polytope and to characterize the class of entanglement witnesses (observables) that 
detect entangled states on opposite sides of the set of separable states. All this serves to motivate a 
classical algorithm which, given the expected values of a subset of an orthogonal basis of observables 
of an otherwise unknown quantum state, searches for an entanglement witness in the span of the 
subset of observables. The idea of such an algorithm, which is an efficient reduction of the quantum 
separability problem to a global optimization problem, was introduced in PRA 70 060303(R), where 
it was shown to be an improvement on the naive approach for the quantum separability problem 
(exhaustive search for a decomposition of the given state into a convex combination of separable 
states). The last section of the paper discusses in more generality such algorithms, which, in our 
case, assume a subroutine that computes the global maximum of a real function of several variables. 
Despite this, we anticipate that such algorithms will perform sufficiently well on small instances that 
they will render a feasible test for separability in some cases of interest (e.g. in 3-by-3 dimensional 
systems) . 



I. INTRODUCTION 

Deciding whether a quantum state, be it physical or 
theoretical, is separable (as opposed to entangled) is a 
problem of fundamental importance in the field of quan- 
tum information processing and is a computationally in- 
tractable problem 0. 

One way to decide that a state is entangled is to use an 
entanglement witness (EW) 0, Q . Much work has been 
done on entanglement witnesses (EWs) and their utility 
in investigating the separability of quantum states, e.g. 
0, 13 ■ EWs have been found to be particularly useful 
for experimentally detecting the entanglement of states 
of the particular form p\ip){^\ -I- (1 — p)a, where \'4)) is 
an entangled state and cr is a mixed state close to the 
maximally mixed state and < p < 1 |^ ■ 

We will show that the set of separable states is not 
a polytope and thus there is no finite set of EWs that 
can detect every entangled state. This work focusses on 
the principle of invoking a sufficiently large subset of all 
possible EWs given the expected values of a set of observ- 
ables. Section HTl summarizes some geometric aspects of 
the set of Hermitian operators and section llTTI reviews the 
geometry of separable states and entanglement witnesses. 
The simplest case of a set of expected values giving rise 
to more than one EW is characterized in section Hvl and 
illustrated by visiting the problem of deciding whether 
a noisy Bell state is entangled. In sectional we apply 
the above principle to the problem of detecting the en- 
tanglement of a completely unknown quantum state and, 
in section IVTI we outline a class of (classical) algorithms 
that search for an EW that detects the state or conclude 
that so such EW exists (given the currently available in- 
formation about the state). 



II. GEOMETRY OF VECTOR SPACE OF 
HERMITIAN OPERATORS 

Let Ti-M^N denote the set of all Hermitian operators map- 
ping C*^(8)C^ to itself. This vector space is endowed with 
the Hilbert-Schmidt inner product {X^Y) = ir{AB), 
which induces the corresponding norm ||X|| = ■\/tr(A"2) 
and distance measure \\X — By fixing an orthog- 
onal Hermitian basis for Ti.M,N, the elements of 'Hm,n 
are in one-to-one correspondence with the elements of 
the real Euclidean space M*^ ^ . Let B = {Xi : i = 
0, 1, ... , M'^N'^ — 1} be an orthonormal, Hermitian basis 



for 



11 A/, TV 



where Xq 



fl. For concreteness, we can 



assume that the elements of B are tensor-products of the 
(suitably normalized) canonical generators of SU(M) and 
SU(N), given e.g. in [|. Note tr(X,) = for all i > 0. 

Define v : Um.n ^ M^^'^'-^ 



as 



viA) 



tr{XiA) 
triX^A) 

tr(XM27V2-l^) 



(1) 



Via the mapping v, the set of separable states Sm,n 
can be viewed as a full-dimensional convex subset of 



{v{a) G 



(2) 



which properly contains the origin v{Im,n) = G 
]gA/ N -1 (pgcall that there is a ball of separable states 
of nonzero radius centred at the maximally mixed state 
Im,n i). 
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Most of the definitions in the rest of this section may 
be fomid in If A e 'Hm,n and A ^ Q and a S M, 

then {x € Ti.M,N ■ tr(Ax) < a} is called the halfspace 
HA,a- The bomidary {x £ Hm.n ■ tr(^a;) = a} of 
i?A,a is the hyperplane TTA,a with normal A. Call two 
hyperplanes parallel if they share the same normal. Let 
denote the interior HA,a \ T^A.a of HA,a- Note that 
_ is just the complement of HA,a- For example, 



H° 



the density operators of an M by A'^ quantum system lie 
on the hyperplane tt/ i, where / is the identity operator. 
Let T>M,N = {p G Hm.n ■ P > 0} n TT/^i denote the 
density operators. 

The intersection of finitely many halfspaces is called a 
polyhedron. Every polyhedron is a convex set. Let D be 
a polyhedron. A set C D is a face of D if there exists a 
halfspace HA,a containing D such that F = D D iTA.a- If 
u is a point in D such that the set {v} is a face of D, then 
u is a vertex of D. A facet of Z? is a nonempty face of D 
having dimension one less than the dimension of D. A 
polyhedron that is contained in a hyperball {x G 'Hm,n ■ 
tr(a;^) = R^} of finite radius i? is a polytope. 



III. SEPARABLE STATES AND 
ENTANGLEMENT WITNESSES 

The set of bipartite separable quantum states Sm,n in 
'Hm,n is defined as the convex hull of the separable pure 
states {|Q;)(a| (8) |/3)(/3| € 'HM,N}a,0, where \a) is a unit 
vector in and |/3) is a unit vector in C^. Let £m,n ~ 
T^M,N \ Sm,n be the set of entangled states. For each 
entangled state p there exists a halfspace Ha. a whose 
interior H'^ ^ contains p but contains no member of Sm,n 
0. Call A e Hm.jv an entanglement witness Q if for 
some a e K 



Sm,n n Ha , 



and £m,n H ^ ^ 



(3) 



Entanglement witnesses A with a = in (j2J) correspond 
to the conventional definition of "entanglement witness" 
found in the literature, e.g. @. 

Entanglement witnesses can be used to determine that 
a physical quantum state is entangled. Suppose A is an 
EW as in Q and that a state p that is produced in the 
lab is not known to be separable. If sufficiently many 
copies of p may be produced, then repeatedly measuring 
the observable A of p gives a good estimate of the expected 
value of A 

{A)p tr(Ap) 

which, if less than a, indicates that p G H"^ ^ and hence 
that p is entangled. Otherwise, if {A)p > a, then p may 
be entangled or separable. The best value of a to use 
in (121) is a* = min|^)(^|g5jj^{(i/)| A j?/))} since, with this 
value of a, the hyperplane tta.u is tangent to Sm.n and 
thus the volume of entangled states that can be detected 
by measuring observable A is maximized. With this in 



mind, define 

a*(A):= nain 

\if>){if>\eSM,N 

if A is an EW. 

Detection of the entanglement of reproducible physical 
states in the lab would be straightforward if there were 
a relatively small number K of EWs Ai such that £m,n 
is contained in 

K 



i=l 



where Oi := a*{Ai). This would imply that Sm,n is 



K 

n 



A,, 



that is, that Sm,n is the intersection of finitely many half- 
spaces. Invoking the isomorphism between 'Hm,n and 
jjA/^w^^ this says that Sm,n is a polytope in R^^'^'-i. 
Minkowski's theorem says that every polytope in R" 
is the convex hull of its finitely many vertices (extreme 
points). Recall that an extreme point of a convex set is 
one that cannot be written as a nontrivial convex combi- 
nation of other elements of the set. To show that Sm,n 
is not a polytope, it suffices to show that it has infinitely 
many extreme points. The extreme points of Sm,n are 
precisely the product states, as we now remind ourselves 
(see also ^3): A mixed state is not extreme, by defini- 
tion. Conversely, we have that 



implies 



which implies that | ('i/'ilV') I — 1 for all i; thus, a pure 
state is extreme. Since Sm,n has infinitely many pure 
product states, we have the following fact, which settles 
a problem mentioned in |l2l | . 



Fact 1. Sm n is not a polytope in 



IV. AMBIDEXTROUS ENTANGLEMENT 
WITNESSES 

Suppose that A is not an entanglement witness but that 
—A is. In this case, an estimate of tr{Ap) is just as 
useful in testing whether p is entangled. We extend the 
definition of "entanglement witness" to reflect this fact: 
Call A € Hm.at a left (entanglement) witness if Q holds 
for some a € M, and a right (entanglement) witness if 



Sm.n n H'LA.-b 



and £A/.iv n ^ _f, 7^ (6) 
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for some 6 G M. As well, for A a right witness, define 

max {(V'lAI^)}. 

\if>){if>\eSM,N 

Note that A is a left witness if and only if —A is a right 
witness. 

The operator A € Hj\/,Af defines the family {TrA.a}ai£R 



of parallel hyperplanes in 



Consider the hyper- 



plane TTA 



which cuts through Sm,n at the 



maximally mixed state Imn- When can tt^ be shifted 
parallel to its normal so that it separates Sm,n from some 
entangled states? If A is both a left and right witness, 
then -KA can be shifted cither in the positive or negative 
directions of the normal. In this case, the two parallel 
hyperplanes T^A,a*(A) ^'iid iTA^'iA) sandwich Sm^n with 
some entangled states outside of the sandwich, which wc 
wiU denote by W{A) := H_A-a'(A) n F_a,-6*(A)- 

Definition 1 (Ambidextrous entanglement wit- 
ness). An operator A € 'Rm,n is an ambidextrous (en- 
tanglement) witness if it is both a left witness and a right 
witness. 

If A is an ambidextrous witness, then p is entangled 
if {A)p < a* (A) or if {A)p > b*{A). We can further de- 
fine a left-handed witness to be an entanglement witness 
that is left but not right. Say that two entangled states 
pi and p2 are on opposite sides of Sm.n if there does 
not exist a halfspacc i?A,a such that i?^ ^ contains pi 
and p2 but contains no separable states. Ambidextrous 
witnesses have the potential advantage over conventional 
(left-handed) entanglement witnesses that they can de- 
tect entangled states on opposite sides of 5a/, jv with the 
same physical measurement. 

Entanglement witnesses can be simply characterized 
by their spectral decomposition. In the following, sup- 
pose A £ HIa/.jv bas spectral decomposition A = 



J2i=o Ai|Ai)(Ai| with Ao < Ai < . . . < Xmn-i- 

Fact 2. The operator A is a left witness if and 
only if there exists k G [0,1, . . . , M N — 2] such that 
span({|Ao) , |Ai) , . . . , |Afc)}) contains no separable pure 
states and Xk+i > Xk- 

Proof. Suppose first that there exists no such k. Then 
|Ao) is, without loss of generality, a separable pure state 
(because the eigenspace corresponding to Aq must con- 
tain a product state), so A cannot be a left witness. To 
prove the converse, suppose that such a k does exist and 
that Afe+i > Afc. Define the real function /(cr) := tr(Acr) 
on Sm,n- Since span({|Ao) , |Ai) , . . . , |Afe)}) contains no 
separable states and Afc+i > Xk, the function satisfies 
/(cr) > Aq. Since the set of separable states is compact, 
there exists a separable state a' that minimizes /(cr). 
Thus, setting a := /(cr') gives Sm,n H Ha a = 0- As 
well, Emm^Ha a 7^ since tr(A|Ao) (Ao|) = Aq < a, and 
so A is a left witness. □ 

Theorem 3. The operator A is a left or right en- 
tanglement witness if and only if (i) there exists k G 



[0, 1, ... , MN - 2] such that span{\Xo) , |Ai) , . . . , |Afc)} 
contains no separable pure states and A^+i > Aj,, or (ii) 
there exists I €z [1,2,..., MN — 1] such that 
span{\Xi) ,\Xi+i) , . . . ,\Xmn-i)} contains no separable 
pure states and Xi > A;_i. 

Theorem 121 immediately gives a method for identifying 
and constructing entanglement witnesses. 

Definition 2 (P arti al Product Basis, Unextendible 
Product Basis |13|^. A partial product basis of C^^ (8) 
is a set S of mutually orthonormal pure product 
states spanning a proper subspace of C^^ (X) C^. An un- 
extendible product basis of C^^ (E) is a partial prod- 
uct basis S of C^^ (E) whose complementary subspace 
(spanS*)^ contains no product state. 

We can use unextendible product bases to construct am- 
bidextrous witnesses. Suppose B is an unextendible 
product basis of C^^ (E) C^, and let B' be disjoint from 
B such that BUB' is an orthonormal basis of E) C^. 
One possibility is the left witness defined by A' as 



\\)eB' 



(7) 



As well, we could split B' into and and define an 
ambidextrous witness A" as 

^" = - E \^l){Xl\+ J2 (8) 

Another thing to realize is that spanS may contain an 
entangled pure state, which can be pulled out and put 
into a (-l-l)-eigenvalue eigenspace of A' . Depending on 
B (and the dimensions M, A''), there may be several mu- 
tually orthogonal pure entangled states in spanS whose 
span contains no product state; let i?" be a set of such 
pure states. Define the ambidextrous witness as 



\\)eB' \\)eB" 



(9) 



This suggests the following problem, related to the com- 
binatorial |l4l | problem of finding unextendible product 
bases: Given M and N, find all orthonormal bases B for 
C*-^ ® such that 

• i? is the disjoint union of Al, B, Ah, 

• spanA/^ and spanA/? contain no product state, 

• span(A/^ U A/j) contains a product state, and 

• minllAil, |Aij|} is maximal. 

Such bases may give "optimal" ambidextrous witnesses, 
which detect the largest volume of entangled states on 
opposite sides of Sm.n- 

The functions a* and b* arc difficult to compute |29l| . 
Thus a criticism of constructing witnesses via the spec- 
tral decomposition is that even if you can construct the 
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corresponding physical observables, you still have to per- 
form a difficult computation to make them useful. How- 
ever, most experimental applications of entanglement 
witnesses are in very low dimensions, where computing 
a* and b* deterministically is not a problem - it may even 
be done analytically, as in the example at the end of this 
section. 

Ambidextrous witnesses represent the simplest case of 
the principle of invoking as many (left) entanglement wit- 
nesses as possible given the expected values of each ele- 
ment of a set X of linearly independent observables, that 
is, the case |X| = 1. In the Sectional we see how this 
principle generalizes to \X\ > 1. 

A simple illustration of how AEWs may be used in- 
volves detecting and distinguishing noisy Bell states. De- 
fine the four Bell states in C^: 

:= (|00)±|11))/V2 
:= (|01)±|10))/V2. 

It is straightforward to show that the Bell states are, 
pairwise, on opposite sides of 52.2. Suppose a left entan- 
glement witness W, with a*{W) = 0, detects jV'^) and 
10^). Without loss of generality, W can be written in the 
BeU basis {\ip+) , !</>+),...} as 



W 



— ei a + bi x x' 

a — hi —£2 X X 

X XXX 

X XXX 



(10) 



for ei and £2 both positive, 
-^d?/'^) ± \4>^)) are separable. Requiring (s+| W |s+) > 
gives 2a > €1 + €2 and requiring (s^| W \s^) > gives 
2a < — ei — £2, which; together, give a contradiction. 
Similar arguments hold for the other pairs of Bell states. 
Define the operators 



A. 
A, 



Both A^ and A^f, are easily seen to be AEWs. It is also 
straightforward to compute the values 



and 



a*{A^) = a*{A^) = -1/2 



b*{A^)=b*{A^)^ +1/2. 



Suppose that there is a source that repeatedly emits 
the same noisy Bell state p and that we want to decide 
whether p is entangled. Define the Pauli operators: 



CTi 

CT2 

0-3 



1 



|1)(0|) 
|1)(0|) 
(|0)(0|-|1)(1|) 



(|0)(0| 

-;^(|0)(1| 
72 



where {|0) , |1)} is the standard orthonormal basis for C^. 
Noting that 

A^ — (Ti ® (Tl — 0-2 (72 
A^ = CTi ® (Tl 4- 0-2 ® 0-2, 

measuring the expected value of the two observables 
(Ti ® CTi and (72 (K) (72 may be sufficient to decide that 
p is entangled because p G £2,2 if one of the following 
four inequalities is true: 



{ai (g) ai) p ± {a2 (E) (T2) p > 1/2 
(cTi «) cri)p ± ((T2 «) cr2)p < -1/2. 



(11) 



If the noise is known to be of a particular form, then 
we can also determine which noisy Bell state was being 
produced. Let \B) be a Bell state. Suppose p is known to 
be of the form p\B) {B\ + {1 ~ p)(T for some a inside both 
sandwiches W{Ati,) and W{A^). With a so defined, one 
of the four inequalities 111|) holds only if exactly one of 
them holds, so that \B) is determined by which inequality 
is satisfied. We remark that, if a and \B) are known, 
knowledge of the expected value of any single observable 
A may allow one to compute p and hence an upper bound 
on the I2 distance between p and the maximally mixed 
state 1/4. This distance may be enough information to 
conclude that p is separable by checking if p is inside the 
largest separable ball centred at //4 



V. DETECTING ENTANGLEMENT OF AN 
UNKNOWN STATE USING PARTIAL 
INFORMATION 

We now consider the task of trying to decide whether 
a completely unknown physical state p, of which many 
copies are available, is entangled. For simplicity, we re- 
strict to p S 1^2,2 but the discussion can be applied to a 
bipartite system of any dimension, replacing Pauli oper- 
ators with canonical generators of SU(M) and SU(N) (or 
any linearly independent Hcrmitian product basis). For 
such p, this problem has already been addressed in |l5j . 
where the so-called "structural physical approximation 
of an unphysical map" [T^ was used to implement the 
Peres-Horodecki positive partial transpose (PPT) test 
0, 01 . While the structural physical approximation is 
experimentally viable in principle, it is currently very 
difficult to do so. Thus, the easiest way to test for en- 
tanglement at present is to perform "state tomography" 
in order to get good estimates of 15 real parameters that 
define p, then reconstruct the density matrix for p and 
carry out the PPT test t^iis matrix. 

An experimentalist has many choices of which 15 pa- 
rameters to estimate: the expectations of any 15 linearly 
independent observables qualify, as do the probability 
distributions of an y 5 mutually unbiased (four-outcome) 
measurements Whatever 15 parameters are cho- 

sen, we assume that the basic tool of the experimentalist 
is the ability to perform local two-outcome measurements 
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on each qubit, e.g. measuring ai on the first qubit and 
a2 on the second. Under this assumption, the scenario 
where the two qubits of p are far apart is easily handled if 
classical communication is allowed between the two labs. 
We further assume, for simplicity, that the set of these 
local two-outcome measurements is the set of Pauli oper- 
ators {(7i}i=o^i.2,3 (defined on page 2)). If di is measured 
on the first qubit and (Tj on the second, repeating this 
procedure on many copies of p gives good estimations of 
the three expectations (ui (81 CTq), (cto ® '^j)^ and {oi ® <Tj) 
(where the subscript "p" is omitted for readability). Let 
us call this procedure measuring aiaj. 

Suppose the experimentalist sets out to solve our prob- 
lem and begins the data collection by measuring cricri 
and then a2(T2- Even though only 6 of the 15 indepen- 
dent parameters defining p have been found, the example 
in the previous section shows that p is entangled if one 
of the four inequalities is true. It is straightforward 
to show that if none of these inequalities is true, then no 
entanglement witness in the span of {cti (g) ai,a2 <8i 172} 
can detect p if it is entangled (s^. However, there may 
be an entanglement witness in the span of 



{(To ® O"!, CTQ ® Cr2, CTl ® CTl, (72 ® a2,(Jl ® CTq, O"! (g) (Tq} 



that does detect p. 

More generally, at any stage of the data-gathering pro- 
cess, if we have the set of expectations {{(Ti(g)aj) : e 
T}, then p is entangled if there is an entanglement wit- 
ness in the span of {ai (g) aj : G T} that detects p 
(TC {{k,l) : fe,/e {0, 1,2,3}}\(0,0)). If the experimen- 
talist has access to a computer program that can quickly 
discover such an entanglement witness (if it exists) , then 
the data-gathering process can be terminated early and 
no more qubits have to be used to decide that p is entan- 
gled. The algorithms described in the next section are 
just such programs. To see this, note that the projection 
^2,2 of ^2,2 onto span{tTi (g) aj : G T} is a full- 
dimensional convex subset of R'^', and the projection p 
of p onto span{CTi g) (Tj : G T} is a point in RI"^' 
such that p ^ ^2,2 if and only if there is an entanglement 
witness in the span of {ciigaj : (i,j) £ T} that detects p. 
Since the following algorithms can be applied to any full- 
dimensional convex set (satisfying certain conditions), we 
can apply them to 52,2- 

We view any such algorithm as an extra tool that an 
experimentalist can use to facilitate entanglement detec- 
tion and minimize the number of copies of p that must 
be measured - essentially, trading classical resources for 
quantum resources. As we saw in the case of constructing 
ambidextrous witnesses, the primary classical resource 
required to invoke a sufficiently large subset of all such 
entanglement witnesses is a subroutine for computing the 
function b* (cquivalcntly, a*). 



VI. ALGORITHMS FOR FINDING 
ENTANGLEMENT WITNESSES BASED ON 
GLOBAL OPTIMIZATION 

Assume that p € 'Di\i^n is a state whose separability is 
unknown. We can handle two scenarios ~ one experimen- 
tal, as described above, and the other theoretical. In the 
theoretical scenarios, we assume that we know the den- 
sity matrix for p; this corresponds to having gathered all 
M^TV^ — 1 independent expected values in the experimen- 
tal scenario. Since the algorithms find an entanglement 
witness when p G £m,n, they could also be applied when 
p is known to be entangled but an entanglement witness 
for p is desired (though one may want to apply the entan- 
glement witness optimization procedure [j| to the result 
of the algorithm, as these algorithms do not necessarily 
output optimal entanglement witnesses). 

Let j be the number of nontrivial expected values of 
p that arc known, 2 < j < M^N'^ — 1; that is, (without 
loss) assume we know the expected values of the elements 
of B' = {Xi, X2, ■ ■ ■ , Xj}. The algorithms either find 
an entanglement witness in span(S') for p, or conclude 
that no such witness exists. For any Y G Ha/ jv with 

Y = X]i=o Di^ii Ifit Y be the j-dimcnsional vector 
of the real numbers yi for i = 1, 2, . . . , j. Define 

Sm,n = {ct : cr G Sm,n}- (12) 

Note Sm,n is a full-dimensional convex set in W ^ prop- 
erly containing the origin (since Im^n is the zero-vector 
inRj). 

Let if be a full-dimensional convex subset of E" which 
contains a ball of finite nonzero radius centred at the ori- 
gin and is contained in a ball of finite radius R. The al- 
gorithms are general and can be used to decide whether 
a hyperplane exists which separates a given point p from 
any given K satisfying these properties. Thus, the clear- 
est way to describe how the algorithms work is to use 
the application-neutral notation of convex analysis. For 
X G R" and 5 > 0, let B{x,5) := {y G M" : ||a;-y|| < 5). 
For a convex subset K C K", let ^(A', 5) := \Jx^kB{x, 5) 
and S{K^ —6) := {x : B (x. S) C A'}. Define the following 
convex body problems |20j |: 

Definition 3 (Weak separation problem for K 
(WSEP(Ar))). Given a rational vector p G and ra- 
tional 6 > 0, cither 

• assert p G S{K,S), or 

• find a rational vector c G M" with||c||oo = 1 such 
that c^x < c^p for every x & K |31j . 

Definition 4 (Weak optimization problem for K 
(WOPT(A'))). Given a rational vector c G K" and ra- 
tional e > 0, either 

• find a rational vector y G such that y G S{K, e) 
and cT^ X < <Fy -\- e for every a; G AT; or 

• assert that S{K, —e) is empty [3^ . 



6 



By taking S and e to be zero, we implicitly define the 
corresponding strong separation (SSEP) and strong op- 
timization (SOPT) problems. Note that by taking K 
to be Sm.n and p to be p (for some state p G 'Dm.n), 
SSEP(iir) corresponds to the problem of finding an en- 
tanglement witness for p (or deciding that one docs not 
exist in the span of B') |33l |: and by further taking c to 
be A, for some A £ Mm.n, SOPT(A') corresponds to the 
problem of computing b* (A) (actually something at least 
as hard, since SOPT(ii') asks for the maximizer of c^x, 
over X £ K, rather than just the maximum). 

We describe oracle- polynomial-time algorithms for 
WSEP(i^), with respect to an oracle for WOVT{K); 
that is, assuming each call to the oracle is assigned unit 
complexity cost, the algorithms have running time in 
0(poly(n,log(i?/J))). We will use "O" to denote ora- 
cles (black-boxed subroutines) for problems, indicating 
which problem via a subscript, e.g. C'sopt(A')- In what 
follows, so as not to obfuscate the main idea of the al- 
gorithms, we ignore the weakness of the separation and 
optimization problems; that is, we assume we are solving 
SS^V{K) with an oracle for SOPT(ii:). 

There are at least two ways to reduce SSEP(j<' ) t o 
SOPT(K). The first method was covered in |2ll [2^: 
the second method, which we give below, is well known 
and may be found in the synthesis of Lemma 4.4.2 and 
Theorem 4.2.2 in |23. For y G R" and e M, define the 
hyperplane -Ky^b = {x S R" : y'^x = b}. 

Definition 5 (Polar of K). The polar K* of a full- 
dimensional convex set K C M" that contains the origin 
is defined as 

K* := {c e K" : c^x < 1 V.t e K}. (13) 

If c S K* , then the plane tTc.i separates p G M" from K 
when c^p > 1. 

Definition 6 (Feasibility problem for K' 
(FEAS(A''))). Given a convex set K' C R", ei- 
ther 

• find a point k' £ K' , or 

• assert that K' is empty. 

Thus, the separation problem for p is equivalent to the 
feasibility problem for Qp, defined as 

Qp -.^ K* n {c : p^c > 1}. (14) 

As outlined in the next section, to solve the feasibility 
problem for any if, it suffices to have a separation 
routine for K' . Because we can easily build a separation 
routine OssEP(Qp) for Qp out of C'sSEP(if*), it suffices 
to have a separation routine CggEp(if*) for K* in order 
to solve the feasibility problem for Qp Building 
CsSEP(Qp) out of OssEP(x*) is done as follows: 

Routine CssEP(Qp)(y): 
CASE: p"^?/ < 1 



RETURN -p 
ELSE: p'^y > 1 

CALL OssEP(K*)(y) 

CASE: OssEP{K'-){y) returns separating vector q 

RETURN q 
else: OssEP(x*)(y) asserts y £ K* 

RETURN "t/ £ Q" 

It remains to show that the optimization routine 
OgOPT(A:) for K gives a separation routine C'gsEp(K*) 
for K*. Suppose y is given to Oq,opt{k)^ which returns 
k £ K such that y'^x < y^k =: b for all x £ K. If 
6 < 1, then C'gsEP(i<-*) may assert y £ K*. Otherwise, 
C'ggEP(K*) may return k, because tt^.i (and hence Tr^.f,) 
separates y from K*: since k^y = 6 > 1, it suffices to 
note that k^c = c^k < 1 for all c £ K* by the definition 
of K* and the fact that k £ K . 

The plane tt^^i is called a cutting plane, and, to solve 
FEAS(A'') with C'ssEP(_fs'')' ^ cutting-plane algo- 

rithm. All such algorithms have the same basic structure: 

1. Define a (possibly very large) regular bounded con- 
vex set Pq which is guaranteed to contain K' , such 
that, for some reasonable definition of "centre", 
the centre of Pq is easily computed. The set 
Po is called an outer approximation to K' . Com- 
mon choices for Pq are the origin-centred hyperbox, 
{x £ M" : -2^ < x^ <2^, I < i < n} and the 
origin-centred hyperball, {x : x'^x < 2^} (where 
2^ is a trivially large bound). 

2. Give the centre ui of the current outer approxima- 
tion P to OsSEP(if')- 

3. If OssEP(if') asserts £ K"\ then HALT. 

4. Otherwise, say C'ssep(a:') returns the cutting plane 
Tic,b such that K' C {x : (Fx < b}. Update (shrink) 
the outer approximation P := P H {x : c'^x < 
for some b' > b; the idea is that the new P has 
about half the volume of the old P (i.e. usually Tr^fc 
passes through lo, or nearby it). Possibly perform 
other computations to further update P. Check 
stopping conditions; if they are met, then HALT. 
Otherwise, go to step 2. 

The difficulty with such algorithms is knowing when to 
halt in step 4. Generally, the stopping conditions are 
related to the size of the current outer approximation. 
Because it is always an approximate (weak) feasibility 
problem that is solved, the associated accuracy parame- 
ter 5' can be exploited to get a "lower bound" V on the 
"size" of K' , with the understanding that if K' is smaller 
than this bound, then the algorithm can correctly assert 
that S{K',—S') is empty. Thus the algorithm stops in 
step 4 when the current outer approximation is smaller 
than V. 

Using the above reduction from SSEP(/v) to 
SOPT(A"), there are a number of polynomial-time convex 
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feasibility algorithms that can be applied (see [23 for a 
discussion of all of them). The three most important are 
the ellipsoid method, the volumetric- centre method, and 
the analytic- centre method. The latter are more efficient 
than the ellipsoid algorithm and are very similar to each 
other in complexity and precision requirements, with the 
analytic-centre cutting-plane (ACCP) algorithm in 
having some supposed practical advantages. 

We refer to |23| and (and references therein) for 
a discussion of details, including computer precision re- 
quirements, for either of the ACCP algorithms arising 
from either the reduction in plj or the well-known one 
given here. 

VII. CLOSING REMARKS 

In terms of attempting to find a practical algorithm for 
the quantum separability problem, the skeptic notices 
that such algorithms appear not to offer any advantage 
over other approaches: instead of having to solve one 
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